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1. Introduction 



We are faced with a paucity of techniques for understanding, quantitatively, the dy- 
namics of strongly interacting field theories. One of the most intriguing, and probably 
the least well understood, proposals has been to study the theory in the so called large 
limit. Here N refers to the size of the gauge group which may be, say, SU(N) or SO(N). 
It is well known that these theories do simplify in the large N limit. The free energy and 
correlation functions have a well defined expansion in powers of jj. Moreover the leading 
term — the planar limit — is thought to capture the essential physics. 

One of the intriguing features of this limit is the existence of a master field |jl|]. This 
notion arises from the observation that the large limit is in some sense a classical limit 
with playing the role of h. The factorisation of products of gauge invariant observables, 
with corrections of 0(1/N) bears this out [§. The master field refers to this "classical" 
configuration that dominates the path integral. Its knowledge would enable one to compute 
gauge invariant quantities without performing the functional integral. We simply evaluate 
them at this point in field space. In fact, it can be argued that by a suitable gauge 
transformation the master gauge field A^{x) can be made space time independent. Thus, 
for example, in QCD4 we need obtain just four "cxo x 00" matrices 

Tantalising as this prospect may seem, it has proved rather elusive to work with. We 
need to understand what these "00 x 00" matrices really are. And then again we need 
effective ways to compute them. 

Recently, it has been possible to get a better hold on this concept @, [H- The 
mathematical formalism of non-commutative probability theory seems to be the right 
setting for the master field. Let us briefiy summarise the relevant facts. 

The master field can be thought of as an operator living in a Fock space generated 
by creation operators obeying no relations. In other words, a Fock space spanned by the 
states 

(1-1) 

where 

ai\n) = 0, a,a] = d,j. (1.2) 

Here the subscripts i can take either discrete or continuous values. In an n-matrix model 
i will run from 1 to 71. A general operator, the master field included, will be built out of 
these di's and a]'s. 
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In the theory of non-commutative probabihty, a very special place is occupied by 
the so called free random variables. They are akin to independent random variables in 
usual probability theory. They are best thought of as the large limit of matrices with 
independent distributions. Thus for an independent n-matrix model the master fields Mi 
corresponding to the matrices Mi are free random variables. They turn out to have a 
simple realisation on this fock space. 

oo 

m = a, + Y.^nlia^i^- (1-3) 

An arbitrary invariant correlation function of the theory is then computable using 

(Tr [MJ'^^M"/ . . . M"^--]) = (^0|M"^M"/ . . . M^J- . (1.4) 

It turns out that the Mn*'* have a rather simple physical interpretation. They are deter- 
mined solely by the distribution for Mi and are the connected n point Green's functions. 
In other words, dropping the subscript i, 

M^ = (tr[M^)eonn.. (1-5) 

Or equivalently the generating function for connected Green's functions 

zM(z) = 1 + ^ M„z" (1.6) 

is the inverse function of the resolvent 

oo 
n=0 

A generalisation of eqn.( |1.3|) to an arbitrary matrix model also exists. We refer the reader 
to P] for details on this as well as equations of motion, examples and more. 

However, it has been difficult to compute the master field in this formalism other than 
in cases where the large correlation functions have already been obtained by conventional 
methods such as the saddle point or the Schwinger-Dyson equations. It might then appear 
that our framework merely allows us to rewrite results of the large theory in a compact 
way. This is not true. As we shall see in this paper the example of generalised QCD2 is 
an instance where we can profitably apply ideas from the theory of free random variables 
to explicitly compute many quantities of interest. Due to the nonlinearity of the theory. 
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the usual methods, mentioned earher, do not appear very tractable. The main notion 
that is exploited here is that of a multiplicative free family of random variables — a 
notion that we will elaborate on below. It allows us to construct the Wilson loop average 
for an arbitrary contour from knowledge of certain infinitesimal ones. We also obtain a 
realisation of the generalised Hopf equation as a collective field equation for the theory. 
This will allow us to deduce the collective field theory in the general case of the cylinder. 
Finally, one can construct the master field strength and gauge field. 

Generalised QCD2 , has been argued to be a theory of interest in the search for 
potential string theory descriptions of higher dimensional QCD. It generalises the heat 
kernel action of QCD2. Alternatively, it can be defined as the general perturbation to 
topological YM2: 

Z = J [DA] [DB] exp(- j (tr(zSF) - (1.8) 

While it doesn't have propagating modes, it is nevertheless more complex in behaviour 
than ordinary QCD. For instance, it exhibits the large N phase transition [Rl on Riemann 



surfaces of almost all genera [10|. It therefore seems worthwhile to have at hand, explicit 



results analogous to |Tl[]. From this point of view the generalised Hopf equation that we 
obtain below would have bearing on a world sheet lagrangian description. In Section 2 we 
outline the notion of multiplicatively free families and its appearance in QCD2 where it 
enables us to compute arbitrary loop averages. We proceed in Section 3 to calculate the 
loop averages < tr[f/'^] > for a simple loop which is then sufficient to obtain a general 
one. In Section 4 we show that the behaviour of eigenvalues of loops are governed by the 
generalised Hopf equation. The collective field theory is then reconstructed. Section 5 is 
devoted to the construction of the master gauge field and field strength. Discussions and 
conclusions comprise section 6. An appendix bears the burden of a necessary combinatorial 
calculation. 



2. Loops as Free Random Variables 

We indicate how a knowledge of infinitesimal loops in generalised QCD2 (ordinary 
QCD2 is a particular case) can be used to compute arbitrary loop averages. The basic 
loop entities that we will be working with are the simple loops. (By a simple loop we shall 
henceforth mean a non-self intersecting contour.) Their importance lies in the fact that a 
set of simple loops, based at one point and non-overlapping, correspond to a family of free 
random variables ( Fig. 1). 
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Fig. 1: Three Simple Loops 

To see this, notice that we could work with the heat kernel action for generalised 
QCD2 i, 

z = J Uvul n ^^[^^]; 

L plaquettes ^2 -j^^ 

Zp[Up]=Y,dRXR[Up]e-''^''^^-, 

R 

where A(i?) is an arbitrary function of the Casimirs. Now the self similar nature of the heat 
kernel always allows us to triangulate the plane such that the contours Ci are the borders 
of some of the triangles. Let us denote the U{N) holonomies along Ci by Ui. Then in any 
computation of averages of products of Ui, we can integrate out (only on the plane) all 
other link variables. We are thus left with an equivalent measure Z = j "DUi Z\Ui], 
where the product runs over all the simple loops, Ui. In other words, Ui have independent 
distributions. 

Actually, the UiS are not just free random variables. They comprise what is known 
as a multiplicative free family [^. Briefly what this means is the following: The product of 
two free random variables with distributions, /Ui and 112 is also a free random variable with 
some distribution denoted hj ni® 112- A one parameter family of free random variables, 
such that ® fxt2 = Attit2? "^i^l be called a multiplicative free family. (Or equivalently 
A*si ® A*S2 = A*si+S2: if we redeflne the parameter t ^ s = logt.) 

Here the area plays the role of the parameter s. In other words, for the two simple 
loops Ci and C2 (flg. 1), Uci{^i)Uc2{^2) has the same distribution as C/c^oC2(^i + ^2)- 
This follows from the self reproducing nature of the heat kernel action ( |2.1D together with 
its exponential dependence on the area. A more direct argument is made in 0. 

This allows us to obtain the distribution for a simple loop of flnite area by starting 
from one of inflnitesimal area. The precise way to do this is to use the S-Transform . This is 
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analogous to the Mellin transform of ordinary probabilty theory, i.e. it is multiphcative for 
the product of distributions. Therefore for our multiphcative free family, Sai{z)Sa2{z) = 
SAi+A2i^)- It then follows that Sa{z) is exponential in A. 

The function S{z) for a non-commutative random variable U is constructed as follows 
i: If 

oo 

(/>(i) = 5^(^)|^7"|o)i^ (2.2) 

n=l 

then construct the inverse function x(-2), i.e. (j){x{z)) = z. The S-transform is defined as: 

S{z) = ^^X{z). (2.3) 

In the next section we shall compute Sa{z) for infinitesimal A and then exponentiate it 
to obtain that for finite A. Then we shall take the 'inverse transform' to finally arrive at 
< tr[[/"] >. 

We conclude the section with a sketch of how arbitrary loop averages are computed. 
This requires two ingredients. Firstly, a geometrical decomposition of an arbitrary loop 
into a 'word' built out of simple, non-overlapping loops, based at one point. And secondly, 
that the latter are free random variables. 

The first was illustrated at length in . Here we shall merely give an instance of how 
it works. 
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Fig. 2: A Loop and its Decomposition into Simple Loops 
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With the labeUing in fig. 2 we can write the contour as 

C =(12) (2T) (13) (34) (45) (54) (41) 

= (12)(2T)(14)(43)(31)(13)(34)(41)(13)(34)(45)(54)(41) 

V ' 

= (12) (2T) (14) (43) (31) (13) (34) (45) (54) (41) (14) (45) (54) (41) 

. . V ' l^-4J 

C/i U2 Us 

(13)(34)(45)(54)(41) (14) (45) (54) (41) (14) (45) (54) (41) 

V V V ' 

U2 Ua Ua 

=UiU2U3U2Ul 

The main strategy is to introduce backtracking or 'thin' loops to peel off the simple loops 
Ui corresponding to the various windows. 

Having decomposed the loop into a word, we can exploit the fact that the Ui are free. 
The loop average for (p.4|), for instance, reads as 

W{C) = < tY[UiU2U3U2Ui] > = < tvUi >< tv[U2U3U2Ui] > 

= < trill > (< trt/2 >^< trill > - < trt/2 >^< trt/3 >< trt/| > (2.5) 
+ < trt/| >< trUs >< trUi >). 

Here we have used a recursion relation for free random variables (see [0) to factorise the 
word into 'letters' of the form < tr[t/"] >. Once we evaluate the latter, we are done. 



3. < tr[C7"] > in Generalised QCD2 

As mentioned earlier, we utilise the S- Transform ( p.3|) ,( p^ . This requires us to eval- 
uate < tr[t/"] > for an infinitesimal loop. Again it is the heat kernel definition that is 
useful. By definition 

W^^Y^^r! I)f/xi?(t/)tr[t/"]e-4^(«). (3.1) 

R 

To evaluate the U integral we can use the fermion representation of U(N) [0. The 
characters are now the Slater fermion wavefunctions. In a second quantised language 
tr[C7"^] becomes a fermion bilinear and 

I VUxR{U)tr[U^] = {0\Y,Bi_^B^\R). (3.2) 
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Here \0) is the ground state of the N fermions and not the vacuum. In terms of Young 
tableaux this is the state with no boxes. \R) is the fermionic state labelled by integers 
{ui}. These integers are related to the number of boxes hi in the ith row of the young 
tableau by = (A^ - l)/2 + 1 - z + h^. 

It is now easy to see that the only representations that contribute in ( |3.2| ) have 



hi =n — m 

(3.3) 

h^=l (2<z<m + l) 

with m any integer between and n — 1. These come with a relative weight of ( — 1)"^. 
Therefore 

n-l 

Wr.= Y, (-l)"^dfl^e-^^(^'-). (3.4) 

m=0 

Our task is simplified in that we need only the large limit of ( |3.4| ) and infinitesimal 
area. In other words, the sum 

A A 

PF„(AA)= 5^(-l)-d«^(l-— A(i?^)). (3.5) 

m=0 

Now A(i?) could, in general, take the form 

HR) = a{fc,}C'{fci+2fc2+3fc3+...}(^) (3-6) 

{ki} 

where C^ki+2k2+3k3+...}{R) is a generalised Casimir. In the continuum theory this corre- 
sponds to ( PT8| ) with 

^{B) = Y,aik.}l[^trB^^''- (3-7) 

{ki} i 

For reasons of clarity and computational efficacy we'll henceforth restrict ourselves to the 
sufficiently general case 

^(^) = E]^^^(^) (3.8) 

k 

i.e. the higher Casimirs. In the path integral this translates into 

$(S)=^afctrS^ (3.9) 

k 

The considerations generalise in a straightforward manner, however. 



We can perform the sum (3^) with A(i?j^) given by (|3.8| ). This is undertaken in the 
appendix. The result is 

iw„(A.l) = l-^^a,n(";*;-2) (3.10) 



The S-transform is now readily computed. From (|2.2|) 

71=1 A; ^ ^ 



(3.11) 



According to ( |2.3| )we need the inverse to ^(j). This is accomplished easily since we only 
need it to lowest order in IS. A. 

Xaa{z) = + ^ E + + ^)'"')- (3.12) 

As argued earlier, Sa{z) for finite area is obtained by exponentiating Saa{z). 

Sa{z) =e^^(^) 

Vp(^)=X^a,(l + M(l + ^)'-' ^3^3^ 
= EA;(afc -afc+i)(l + z)'=-i 

A: 

We are now in a position to calculate < tr[t/'^] >, given the explicit form of ( |3.13| ). From 
(p^), we conclude that 

< tr[f/"] >=^ I (t){e-'^)e'''^de. (3.14) 



27r . 

We make the change of variables x(z) = e~*^ exploiting the fact that (/)(x(^)) = 2^- Finally 
after an integration by parts we reach 



:1 + l)-e-"^^fc"^(^+^")(^+")* 

n / 27rz 2; 



(3.15) 



This is the expression that generalises the Laguerre polynomials that appeared in the 
usual QCD2 [jl^. They are, of course, recovered by restricting ourselves to the case with 
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dk = Sk,2- We also note that the ai dependence is only through a multiplicative factor of 

A 

e~""^^. The first few averages are 

< trU >=e"^^fe"^< tr[C/^] >= (1 - A ^ afc(A; - l))e"^^ 

< tr[U'] >=(1 - 3-Y,Mk - m + 2) + 3—[J2ak{k - l)]^)e-'^ 

k k 

The general structure is 



<tr[t/'^] >=P„(A)e"''^Sfc"^ (3.17) 
where PniA) is a polynomial of degree n — 1 determined by 

1 rl^~^ A fe 

P^(vl) = ^— -T-r -[x-e-"^^^'=^"'=-"'=+^)(" I , . (3.18) 

^ ^ (n- l)!da;"-i^ ^ '^=1 ^ ^ 

As argued in the previous section this is all we need to compute the loop average of an 
arbitrary contour. 

4. The Generalised Hopf Equation and Collective Field Theory 

The Hopf equation arises in the collective field theory description in QCD2 fl^l • It 



governs the behaviour of the eigenvalue density of Wilson loops |jT^ . In this case we shall 



find an analogous equation — the generalised Hopf equation |p. We shall demonstrate 
that 

F{e, A) = (t>A{e''^) = e'^R{e'\ A) - 1 (4.1) 

obeys the equation 

OF i ^^^dF , 

with (p{F) as in ( |3.13| ). Here i?(e*^, A) is the resolvent 

00 

R{(^ A) = J2< tr[C/(A)'^] > (-(^+^\ (4.3) 

n=0 

We also note that the eigenvalue density of loops a{0, A) is given by 

ReF{e,A) =Tia{e,A) --. (4.4) 
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If we define 0{z,A) sucfi tliat e^^(^'^) = XAiz), then F{e{z,A),A) = z. Tfiis follows 
from ([4 .11 ) and the definition of Xa{z) as the inverse to (p{j)- Therefore 

dF{9{z,A),A) I aP I dF I dejz^A) ^ 

dA 1^ dA de dA ^ ^ 

Knowing, as we do, Xa{z) for generalised QCD2- we see that 

e{z,A) = '-^^{z)-iH^). (4.6) 

This, together with ( [4 .51) enables us to arrive at ( [4 .21) . 

Once again, the hopf equation of QCD2 is regained in the case ak = Sk,2 after a minor 
shift. Like the Hopf equation, (|4.2| ) is integrable as well. In fact, it is easily checked that 

F(^,A) = Fo(^-^(^(F(^,A))) (4.7) 

is the (implicit) solution of this equation. Here Fo{9) = F{6, A = 0) is the necessary initial 
condition. The solution is also completely determined once we specify, say, a{6, A = Q) 
and cr(6', A = 00). 

One might imagine obtaining this equation from the analogue of the loop equations 
for generalised QCD2. These highly non- linear equations would have to be derived from 
a series of Schwinger-Dyson equations for loops. This does not seem to be a very easy 
approach. 



Alternatively, one might imagine a derivation of (|4.2|) from a collective field theory 
description |T^ of generalised QCD2. Logically, one would start from the fermionic de- 
scription. The hamiltonian is a rather complicated polynomial in the operators ^^^r- One 
would have to then bosonise them in their second quantised form and look at the classical 
equations of motion. Here, we have seen it arise much more simply. 

In fact, though it was derived here for the plane, we can conclude that the equation will 
be true on the sphere, or more generally, the cylinder. This follows since these are local 
evolution equations which are sensitive to the global geometry only via the initial/final 
conditions. Thus, on the cylinder we would specify the boundary holonomies or eigenvalue 
densities a{9,A = 0) and a{6,A = Aq) [Aq is the area of the cylinder). Then a{6,A) 
would be determined from ( [4.2| ) with F{6^ A) a complex function such that 

F{e,A)='Ka{6,A)-- - iv{e) (4.8) 
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This can be interpreted as a general kind of fluid flow with density a and velocity v. 

But now we can turn the logic around and derive the collective field theory Hamil- 
tonian whose classical equation of motion is ( [4.2| ). It is convenient to redefine variables 
to 

G{e, A) = l + F{e, A) = na{e, A) + l- ideU{e, A) (4.9) 
where 11 is the canonical momentum conjugate to a. Then we have the Poisson bracket 

{G{d, A),G{d', A)} = -i2nde6{e - 9'). (4.10) 

The generalised Hopf equation now reads as 

|^ + l5:Ma.-a.+,)G-^^ = 0. (4.11) 

k 

This is the Hamiltonian equation of motion 
dG 
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(4.12) 



This is thus the Collective field hamiltonian which can be expressed in terms of cr{6,A) 
and Il{6,A). Our considerations thus far have been in euclidean space. The minkowski 
version reads as 

Z^k + l J 2 2 ^4_;^3^ 

P± =deU{e,A)±7ra{e,A). 



The usual D = 1 matrix model hamiltonian 15 



[ d9amdell{9f + ^a{ef} (4.14) 



corresponds to taking ak = Sk,2- 
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5. The Master Field 

Given tliat we liave computed arbitrary loop averages, we expect to have enough 
information to construct the master field strength. To actually do this we will first obtain 
master loop operators U{AA) for infinitesimal loops. These take the form 

oo 

k=o (5.1) 



<tr[t/"] >=^0|t/"|0 

where we need Uk only to lowest order in AA. Since U{y) = ^ + Y^'kLo^k+iV^ is the 
inverse of the resolvent, using (4J.) and (|2.2D we have 

R{^) = (1 + z)x{z) ^U{{z + l)x{z)) = (5.2) 
Therefore ( |3.13| ) implies 

U{y)=-{l + z{y)y, y = ze^^^'^ 



y 

oo 
k=0 

This is easily solved to lowest order in A to obtain 



(5.3) 



AA 

z{y)=y{l-^cp{y)). (5.4) 



So the master loop operator U{AA) is 



U{AA) ^d+l-—ip{d^). (5.5) 
2 



Let us now construct the operator H such that = U. This is best done as follows. 

2 

C/" = e*"^ = l + mi?-^^2 + .... (5.6) 

Therfore H is identified as the linear term in n in U"^. Again matters are simplified in 
having to keep only terms of 0{AA). Then from (|5.5| ) 

A A 

AA ^^-^^ 



2 

m=0 
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First look at the term of order A A. Terms in this expansion of the form d^^a"^ do not 
contribute to in leading order in AA. We group the other terms where multiplies 
on the left and look for the contribution linear in n. This gives us the term 



r=0 



T^^f^i^'y (5-8) 



In this expression we can drop the terms purely of the form a"^. Since there are already 
such terms of 0(1) and the 0{AA) only contribute to higher order. The terms can 
be easily read from the A = limit when U = d + 1. Putting it all together 

AA 

H = - i{log{l + a){l + — aV(a^} 



I II' /\ /I 



r 2 ^ " (5.9) 



r=l n 

. ^ (^ + 1) \n + r) 



We notice that cq = 0. 

Matters are simplified if we choose axial gauge Ai — Then the holonomy for an 
infinitesimal loop becomes 

U = Pexp(z ^ A^dx^) = exp{idiAoAA). (5.10) 

We can then identify the large limit of the two matrices: 

FAA = diAoAA = H (5.11) 



with H as in ( |5.9|) . Since cq = 0, F has a vanishing 1-point function which is expected. 
The presence of the AA is a consequence of discretisation of the action. This ensures that 
correlation functions of the field strength F do not have delta functions in them. In this 
gauge the master field Ai = and Aq is determined by (|5.11|) and (|5.9| ). 



6. Discussions and Conclusions 

In this work we have applied the formalism of the master field to explicitly compute 
the large limit of many quantities in generalised QCD2. This was meant to illustrate the 
utility of the framework of non-commutative probability theory in actual computations. 
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As we have seen, it is not easy to see how to obtain the coUective field theory description 
or the expressions for loop averages and gauge fields with the usual means available to us. 
It is worthwhile to make some remarks regarding them here: The collective field equations 
of motion of Section 4 are related to the non-commutative probability distributions in 
the same way as the heat equation is related to the Gaussian distribution. This gives 
us an inkling as to the mathematical relation of the collective field description and the 
underlying matrix model. It might help to understand the relation better so as to arrive at 
tractable collective field descriptions of realistic theories. On a technical level, we note the 
simplification gained in going to an operator construction of the master field. In Section 
5 we were able to go from loop averages to gauge fields, simply by performing standard 
operator manipulations like taking the logarithm. Analogous operations in a large N 
matrix description do not appear easy. 

We have, of course, in this paper, exploited a very special property of this system, 
namely that of being multiplicatively free. This would not straightforwardly generalise to 
higher dimensions. Nevertheless, the lesson that is perhaps to be drawn is that techniques 
can be developed in the operator framework for the master field which can take us beyond 
conventional approaches. 

Acknowledgements: I would like to thank David Gross for suggesting the appli- 
cability of the master field framework in this context, as well as for useful discussions. I 
would also like to acknowledge Ori Ganor for helpful conversation. 

Appendix 

Here we undertake to perform the sum 



m=0 

to leading order in large A^, over representations Rm of ( |3.3| ). The dimension of the 
representation is given by 




n-l 





(6.2) 
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A(i?^) is given by (3^), and the higher Casimirs are in general given by 



N 

1=1 



(6.3) 



In our case ( |3.3| ), this can be computed to be 

m + l7m+l 

= ^[(A^ + n-m - l)''~^{n - m - 1) + {N - m)''"^m] 4^ 



n 



n 



j[{N + n-m-l)'' - {N - mf - N{k ^k- 1)]. 



Using (|6.2| ) it is easy to show that the term in ( |6.1|) is indeed 1. The non-trivial part 
is the 0{/S.A) term which reads 



n-l M n )\m) 

k>0 m=0 ^ / \ / 



(6.5) 



X [{N + n-m~l)'' - {N - m)^ - N{k ^k- 1)]. 



For fixed k consider the (A^ — m)^ term. Expanding in powers of A^, the sum can be 
expressed as the coefficient of y^~^ in 

fc /; \ n— 1 / 1 \ rl 



r=0 ^ m=0 ^ ^ 



Z—— Z „iv-i 



n — 1 \r J dz 

r=0 ^ ' 



z=y,y' 

(6.6) 



Next we write {z-^Y = ^i^i Cfz -j^ , where C[ are coefficients determined by the 
recursion relation 

c[ = ic;-^ + cizl- c; = 1. (6.7) 



Thus (^^) reads as 

-l)\yN-.. (6.8) 



k 

n 



n 



I E(-l)'^ + y)^ Yl CI{-iyy\n - l)(n - 2) . . . (n - 01, 
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Since to leading order in we only need the term that goes as A^'^, it suffices to look for 
the coefficient of A^'' in 

(1 + yfy-l,"-' = + J = (TTT^ + + o(n-^). (6.9) 



Thus only the / = r and I = r — 1 terms contribute in the sum over / in ( |6.8|) . Therefore, 
using = 1 and C^_i = ^-^^^y^, the leading contribution to ( |b.8[ ) simplifies into 



k 



r=0 



It can be checked that the (A^ + n — m — 1)'^ term in ( |6.5| ) gives a contribution of opposite 
sign but same magnitude as ( |6.10| ). Therefore the full contribution from (|6.5|) reads as 

A;>0 ^ ^ ^ ^ A;>0 ^ ^ 

This is the result used in ( |3.10| ). 
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